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Abstract 

We present an approach to higher point loop integrals using Chinese magic in the vir- 
tual loop integration variable. We show, using the five point function in the important 
e+e~ — > // + 7 process for ISR as a pedagogical vehicle, that we get an expression for it 
directly reduced to one scalar 5-point function and 4-, 3-, and 2- point integrals, thereby 
avoiding the computation of the usual three tensor 5-pt Passarino-Veltman reduction. We 
argue that this offers potential for greater numerical stability. 
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With the advent of the LHC, we enter the era of precision QCD, by which we mean 
predictions for QCD processes at the total precision tag of 1% or better. This is anal- 
ogous to the per mille level era of EW corrections at LEP energies. Radiative effects 
at the level of 0{al) have to be controlled on the QCD side and those at the level of 
0{aLas), 0{a^L'^) on the QED(8)QCD and QED sides have to be controlled system- 
atically, both from the physical precision standpoint and from the technical precision 
standpoint, in order to optimize physics discovery at the LHC0. In Ref. [1], we have de- 
veloped a platform for the realization of such corrections ultimately on an event-by-event 
basis based on exact, amplitude-based resummation of QED and QCD together, wherein 
residuals for hard photons and hard gluons are simultaneously calculated order-by-order 
in perturbation theory in powers of a and as- These residuals, which are infrared finite 
and, for hadron-hadron applications, coUinearly finite require then exact evaluation of 
higher point and (higher) loop Feynman diagrams in an appropriate reduction scheme for 
any attendant tensor properties as first developed systematically in Ref. [2], for example. 
Recently, alternative approaches have been developed in Refs. [3lll] to deal with the grow- 
ing complexity of the method in Ref. [2] as the number of legs beyond four and/or loops 
beyond one increases. Here, we focus on higher point one-loop function^. 

It has been demonstrated that n-point functions, for = 1, ■ ■ ■ , 4, at one-loop, re- 
duced to scalar functions using the method of Ref. [2], are tractable for fast MC event 
generator implementation for arbitrary masses and kinematics for high energy scattering 
processes P422]. It has also been demonstrated [231426] that, at one-loop, higher point 
scalar functions can be reduced to sums of four-point scalar functions. In Refs. [M1I271 - I30] 
representations of the scalar four-point function that cover arbitrary masses and the mo- 
menta relevant to most high energy collider applications have been given and these are 
suitable for fast MC implementation. Thus, when one is discussing higher point functions 
at one-loop, we can consider, at least for most collider physics applications, that the 1, 2, 
3 and 4 point functions at one-loop are known in a practical way so that the main issue 
can be considered to be the representation of the higher point functions in terms of these 
known functions. 

When we consider any higher point function, two of the most important aspects of 
any reduction procedure for recasting it in terms of the "known", lower point functions 
are its numerical stability and its usefulness for Monte Carlo event generator realization, 

as we have in mind for our residuals /3n,m in Refs. [1] for example. Given the simplification 
that has been shown for the "Chinese magic" polarization scheme [3TI433] for real emission 
of massless gauge particles in such functions, it is natural to seek further simplification 
and numerical stability in the virtual emission and re-absorption processes as well by 
exploiting the same scheme. It is this that we pursue in what follows. 

For the reader unfamiliar with the "Chinese magic" polarization scheme for massless 

^Here, L denotes the typical big log for the process under discussion. 

^See Refs. [5] for some recent progress on the higher loop functions with an eye toward their use in 
the MC realization of the approach in Ref. [T] 
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gauge bosons, which is historically associated to the preprint in Ref. |3T], the key obser- 
vation is that the gauge invariance of the attendant massless gauge theory allows one to 
use an attendant set of polarization vectors which, when the chiral forms of the respective 
spin I charged particles' wave functions are used, eliminate radiation from one entire side 
of a charged line and, simultaneously, simplify considerably the calculation of the part of 
the amplitude that remains, almost like "magic, hence the name. This is possible because 
of a representation of the respective polarization vector for helicity and 4-momentum 
kj, e'^ , as a matrix element of the Dirac gamma matrix, 7^, between the spinor of helicity 
and four-momentum k^, \k^\^ >, and the massless spinor state < pA^|, = 0, up 
to a normalization factor, so that the Chisholm identity (see eq. fll2p below) reduces the 
Feynman rule factor e^^7^ at the respective interaction vertex to the simple expression 
2[|/c^ — A^ >< p — A^l + |pA^ >< A;^A^|], up to the same normalization factor, which causes 
one side of a line of the real radiation terms to vanish if p is set equal to the external 
4-momentum entering(leaving) that side of the respective line. The remaining terms are 
then expressed in terms of simple spinor products which lend themselves to easy evalua- 
tion piHSB] . This gives a 'magically' shortened expression compared the usual Cartesian 
representation of the polarization vector with the squared amplitude modulus evaluated 
using traces over the fermion lines. We illustrate this below here as well. 



Specifically, we will use the conventions of Ref. [T0J|34] for spinors and polarization 
vectors, which are derived from the work of [STllSSj . The 5-pt function which we want to 
analyze in these conventions as our prototypical example is shown in diagram (c) in Fig. [H 
It has many applications in collider precision physics. When combined with diagrams 
and 1(b) it generates a gauge invariant contribution to the ISR for e~^e~ — ff+'y, f 7^ 
for example, and it is a part of such a contribution to mm — )■ p/2 -|- G (in an appropriate 
color basis), etc. Such applications and their attendant phenomenology will be taken up 
elsewhere. [SB]. Here, we focus on the use of Chinese magic in the loop integral in Fig. 
1(c) to illustrate what simplifications are possible. 

More precisely, by the standard methods, we need the following Feynman integral 
representation of Fig. 1 (c) 

■^A.A.AiA^A, - W Pi P2 l')^ ^^n)^ ^{q + p, - k)^ - mj + ie){iq + pi)^ - mj + ie) 



{{q+Pl+P2 - kY - M2 + +p'J2 _ + ^g)(^2 _ ^2 + 



+ . 



(1) 



where we have defined massless limit coupling factor 

C = C({Ai}, {A;.}) = QieG^G'\v[ + a[\2){vi - ai\i){v'^ + a'^\'^){v2 - aaA'^) (2) 
with the couplings Qie, G and G' for the 7, Vi and V2, respectively. In the usual 



numerical realization of the amplitude in Fig. [T] as it relates to bhabha scattering can be found in 
Ref. 135]. 
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(a) (6) (c) 

Figure 1: ISR 5-point function contributions with fermion and vector boson masses 
ruf, rriB, f = 1,2, B = Vi,V2 and with four momenta Pi, k as shown, with Q = pi + p2. 
Radiation is shown from the initial state hne with electric charge Qie where e is the 
electric charge of the positron - here pi is the incoming fermion 4-momentum, p2 is the 
incoming anti- fermion 4-momentum. When the quantum numbers allow it, the crossed 
graphs for the internal vector boson exchanges must be added to what we show here. 

Glashow-Salam- Weinberg- 't Hooft-Veltman [37] notation, v(a) represents vector(axial- 
vector) coupling. The ellipsis in ([1]) represent the mass corrections needed to correct the 
massless limit used for C{{Xi}, {A^}). They are not necessary to illustrate our method and 
they will be restored elsewhere [36]. To get the loop integral in terms of Chinese magic, 
we take the following kinematics as shown in Fig. [D 

Pi = {.E,pz) 

P2 = {E, -pz) 

— P4 = {E',p' {cos 9[z + sm9[x)) = p[ 

k = {k^, k{cos 9^z + sin 9^(cos(f)^x + siiii(f)yy))) 

-Pa -P3 + k = pi+p2 = {^^s, 0) 

-P3 = P>2 (3) 

with k^ = k, y/s = 2E. Here, we introduce the alternate notations p'l = —p^, p'2 = —ps 
for cosmetic use entirely. We now introduce the two sets of magic polarization vectors 
associated to the two incoming lines: 

^^^^^^ V2u.Ak)uAPy ^ ""^ v^^_.(fcK(C)' ^ ' 

with = and ( defined in Ref. [lOllM], so that all phase information is strictly known 
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in our amplitudes: the two choices for /3 are such that its space-hke components have the 
directions of the two incoming beams in the initial state. We take the the basis of the 
4-dimensional momentum space as follows: 

ii = {E,Ez), h = {E,-Ez) 

V2<i2- > 
= ^^'^''^ (5) 

< ^2 - > 
= ^(^-^^) 

where we use the obvious equivalence \ia >= u{tj^ in the notation of Refs. [3T1431] . The 
important point is that all four of these basis 4- vectors are light-like with if = 0, i = 
1, • • ■ ,4. They therefore can participate in Chinese magic. 

To illustrate explicitly this latter point, consider the definite case Ai, A2, X'l, Ag, A^ = 
+, — , +, — , +, as all other choices for the helicities behave similarly. We write the loop 
momentum as 

q = aiii (6) 

with summation over repeated indices understood. The coefficients ai are readily deter- 
mined as 





qh 




2^2 




qii 


a2 = 


2^2 




qh 


"3 = 


E2 " 




i 


04 = 





(Di-Do-M^J/s 



'2s 

where we define the denominators as 



(S)3-D2-s + 2p2fc + M4)/s 

(7) 

[c,D,- + csM^^ + C6(M4 + 2p2k -s) + C7{2kpi)] 



E2 



2)o = - M^^ + te 

'^1 = iq + Pi)^ - ml + ie 

2)2 = (g + Pi - ky -m\ + ie 

^3 = iq + Pi+P2- kf - M4 + ie 

2)4 = (g — PiY — m\ + ie 
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so that the expansion coefficients {cj} are 

, ,csc6'^e*'^^ CSC 611 e*"^^ esc 61^ cot 6'1 e*"^^ - cot 61^ , 
Co = CSC 57^^ + —7^ -) 



Ci = CSC I 



, CSC 6[e"^^ CSC 9[ cot 9[e"^^ — cot esc 



, — CSC 6', 6*"^^ CSC 6^^ cot 9\ e'^^ — cot 6^^, esc 9^ 

C2 = CSC(j)J — ^ + T . 1 ^ ^ 



C3 = CSC I 



C4 



Cq = CSC I 



, CSC 9[e'"l'^i CSC cot 9'ie"^'^' — cot 6^^ 
CSC 9[e'"^'' 

, CSC 6''i e*"^^ CSC 6*1 e*"^^ esc 9^ cot e*<^^ - cot 6*^ 

, CSC 9[e''^^' CSC 6'^ cot 6';e*'^^ - cot 9^ 



(9) 



CSC 6*^ 

C7 = -CSC0^— 

Thus, the {cj} are determined exphcitly by the cms kinematics that we use. The con- 
sequence to note is that the Chinese magic now carries over to the loop variable via the 
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(10) 



identity 

2 

+ «3 , (14- ><h-\ + >< h + I) 

<P2- \Pl+ > 

+ a4 ^ (1^2+ ><ii + \ + \ei- >< h - 1) 

<P2 + \Pl- > 

2 

= >< Pj + \ + \Pj- >< Pj - 1) 

V^^ /I II IN 

+ "3- ^ — (b2- >< Pi - I + bl+ >< P2 + I) 

<P2- \Pl+ > 

V^^ n II IN 

+ a4 i\p2+ >< Pi + I + bi- >< P2 - I) 

<P2 + \Pl- > 

2 

= Yl >< Pi + I + bi- >< Pi - I) 

i=i 

+ "3(1^2- >< Pi - I + bl+ >< P2 + I) 

+ a4{\P2+ >< Pi + I + |pi- >< P2 - I) 

where we work in the massless hmit for this numerator algebra so that we take £i = 
pi, £2 = p2 in (fTOl) . Here, we defined as well 

. _ V2E as 

= as ; = 1= 

<P2-|pi+> V2 

. _ V2E a^ 

= ^4 ; = — 

<P2 + |pi-> V2 

From the standpoint of efficient and numerically stable MC event generator realization of 
the correction in Fig. [H the explicit form the aj cannot be stressed too much. 



Upon introducing the representation ([10|) into the numerator, N, of the integrand in 
}T]) we get, from the standard identities 

^ [\iiX.f >< kX^\ + \k - >< ii - X^\], 



< k — A-yl^lA-y > 

\ >= 2 

h = ><ii + \ + >< ii 



Y < hXhp\i2X >= 2[|£i - A X £2 - A| + \i2X >< iiX\], ^^^^ 
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the reduction 



N 




{{Ai <P2 + \p'i- >< P2 - \P2+ > +A2 <P2 + \p'l- >< P'2 - \Pl+ >) 



(A3 <P2 + \Pi- >< Pi - \P1+ > +A4 <pi + \p[- >< p[ - \pi+ >) 
+ 64(^1 <P2 + \P1- >< p'2 - b2+ > +^2 < P2 + bl- >< J>2 - bl+ >) 
(A3 < P2 + \Px- >< P2 - \Pl+ > +A4 <Pl + \Pi- >< P2 - \Pl+ >)}, 



(13) 



where we defined 



A 



1 



0(4 < Pi + \k~ > +a2 < P2 + \k- >, 

(1 + tti) < pi + \k- > +^3 < P2+ \k- > 

"2 < Pi - \P2+ >, A4 = a4 < pi - \p2+ > 



A2 
As 



(14) 



for the magic choice /3 = pi. Note that the 'magic' has killed all but one set of the 
terms with three factors of the virtual momentum expansion coefficients and that, in 
the numerator of the propagator (before) after the real emission vertex, it has eliminated 
the terms associated with {^i) /c as well as half of the terms in the respective virtual 
momentum expansion in former case. While we have eliminated a large fraction of the 
possible terms on the RHS of (|T3|) . one can ask how it compares in length with what one 
would get from the usual approaches of taking traces on the fermion lines. To be specific, 
in the traditional method that leads to traces on fermion lines, one needs to compare the 
length of 2^Ai*^Ai^^^^ where is the respective Born amplitude that would interfere 
with the one-loop amplitude to create the one-loop correction to the respective cross 
section. In the Chinese magic representation, we get immediately that only radiation 
from the anti-particle (^2) incoming line contributes with the simple result (repeated 
indices are summed and s' = [pi + p2 — k^ as usual) 



so that computing 23?A^^A^*^^'^^ just involves multiplying in f|T3|) by the complex 
conjugate of this simple expression and taking twice the real part. If we proceed with the 
usual trace on the fermion lines method, one needs the trace of two sets of terms with 10 
Dirac gamma matrices multiplied by a factor with the trace of 6 Dirac gamma matrices: 
this means one has 2-9-7-5-4x5-4 = 2520 x 20 = 50, 400 terms, each of which requires 
Passarino-Veltman reduction of 3, 2, and 1 5-pt tensor integrals. In Ref. [38], another 
approach that leads as well to traces over fermions is used in which one first expands the 
amplitude under study in a gauge invariant tensor basis with scalar coefficients and uses 
Chinese magic-type representations of the helicity states to express the attendant 

helicity amplitudes in terms of these invariant scalar coefficient functions. The key step 
is the use of projection operators, V{X) in the notation of Ref. [38], which project out 
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the scalar coefficient X. To compare with our approach, we observe the following: the 
Born amplitude tensor structure is one of the tensor structures in the respective expansion 
basis and to project its coefficient the respective projection operator evaluates a linear 
combination of the trace on the fermion lines of the hermitian conjugate of this Born level 
tensor structure in product with the Feynman amplitude and the traces on the fermion 
lines of the hermitian conjugates of the other tensor structures in product with the same 
amplitude. Thus, our counting of terms given for the evaluation of 2'SiAi*B-M.^^^^ using 
the traditional traces on fermion lines gives a lower limit to the number of terms that 
would be generated by the methods of Ref. [38] for our calculatioij^. Looked at this way, 
we can appreciate better the great simplification that ( !T3|) represents. It follows that this 
form of in ( lT3l) has efficiently reduced the problem of reduction of the 5-pt function 
with three, two and one tensor indices (index) in the Passarino-Veltman formalism to the 
problem of a single scalar 5-pt function and lower 4, 3 and 2 point functions with the 
coefficients already explicitly expressed in terms of the cms kinematic variables that are 
so crucial to efficient MC event generation. Efficient MC event generator realization of 
the latter functions is known [61422]. where it is understood that one uses the results in 
Ref. [25H2B] to express the scalar 5-pt function in terms of scalar 4-pt functions using our 
explicit kinematics above. These last remarks are made more manifest when one notes 
the introduction of the result for in (fT3|) into the integral in (JT]) leads to the integrals 

all of which are known from the lower point functions we advertised when the results for the 
representation of the scalar 5-point function in terms of 4-point functions in Refs. [21H2B] 
are usecif]. We get a bonus: no evaluation of wave functions at complex momenta is 
required here. What we have done is rigorously a result of Lagrangian quantum field 
theory and it therefore can serve as a cross check on methods that may not obviously so 
be. Evidently, the method we have illustrated can be used for any higher point function. 



■^For example, let us take the example discussed in Ref. [38], using their notation, of 
<]iP2, A2) ^(pi, Ai) — >■ 7(p3, A3) 7(p4, A4), where we focus just on the one-loop correction from the Gross- 
Wilczek-Politzer [55] QCD theory with direct analysis for the respective 4-pt box graph in which a 
gluon is exchanged between the incoming quark((7) anti-quark(g) pair "before" they annihilate to the 

two photons. For the helicities H h-l- for the quark, anti-quark, 7(^3), 7(^4), respectively, the helicity 

amplitude is proportional to the An scalar coefficient in Ref. |38| . Evaluation of the projection operator 
for All on the box graph requires the trace for a product of 12 Dirac gamma matrices, which generates 
11 • 9 • 7 • 5 • 3 = 10, 395 terms, and this has to be done 5 times (there are five scalar coefficients) for a 
total of 51, 975 terms. This is just a 4-pt function. The same calculation using our methods generates a 
formula smaller in length than that in eq. (j20p in the text. 

^ Since the integral in ([l} is manifestly UV finite, we do not need to specify what regularization is 
used for the two point functions because only UV finite combinations of them can occur here while the 
wave functions arc all in 4-dimensional Minkowski space. Note also that the standard trace over fermion 
lines would also lead to results equivalent to that in but as we have seen above it would necessitate 
evaluation and simplification of much longer expressions in general to compute the attendant transition 
rate for the process. 
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At this point, while we have shortened considerably the respective amplitude and 
have removed the Gram determinant type factors in the tensor reductions, we are still 
subject to the Gram determinant-type denominator factors in the results in Refs. 
for the representation of the 5-point scalar function in terms of 4-point scalar functions. 
We have found that these are in general still too numerically unstable for realization in the 
amplitude-based exact resummation MC event generators such as those in Refs. [9]. Thus, 
we replace the representation from Refs. pH - [26] of the needed 5-point scalar function here 
as follows. 



We start from the basic identity 



g2 = Do + m2 - te 



= 2aia2^i^2 + 2^30444 ^^^^ 
s 

= saia2 + -«3a4- 

Dividing by Dq ■ ■ ■ 2)4 and integrating over d^q we arrive at the following representation 
of the required scalar 5-point function (we use the notation of Ref. [26] for Eq itself): 



Eo{pi,p2,p3,p4,mo,mi,m2,m3,fh^) = {- A(0) + ^^Po(13) - Co(12) - Co(03) + Co(02) 

+ (M4 - s + 2p2k){Do{l) - DoiO)) - M2^Po(3) - Do{2 
- l^[Ari,o(I)o(l) - DoiO)) + 2Api,o(I?ii(l)p(l)i 



- L'ii(0)p(0)i + I)i2(l)p(l)2 - A2(0M0)2 + Di3il)m)3 

- 1)13(0)^(0)3) - I)o(l)p(l)4 + I)o(0)p(0)4 

- 2M2^ (Do(1) - ^o(O)) + Ar3,2(I)o(3) - Do(2)) 

+ 2Ap3,2(L'ii(3)p(3)i - Dii (2)^(2)1 + Z)i2(3)p(3)2 - 2^12(2)^(2)2 
+ 1)13(3)^(3)3 - I)i3(2)p(2)3 - I)o(3)p(3)4 + 1)0(2)^(2)4) 

, 4 

+ 2(m4 - s + 2p2k)iDo{3) - Do(2))] - -^[^\cj\\Co{j, j + 1) 

j=0 

+ Arjj+iDoU) + 2APj^j+i{Du{j)p{j)i + Di2U)pUh 

4 

+ Dr,{j)p{j)3 - Doimjh)) + 2(^ ^{c,c*)Co{ij) 

i<j 

4 

+ J^3f?(c,(c^M2^ + c^(m4 - s + 2p2k) + c*,{2kpi)))Do{j))] }/Ceo, 

j=0 

(18) 

where we have the identifications 

Pi = Pi, p2 = Pi-k, p3 = Pi+P2-k, Pi = p'l, rho = My,, rhi = mi, 7712 = mi, 7713 = M^j, 7714 = m2 
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and where the coefficient Ceq is given by 

Ce, = -ie + ^t^M2^(M4 - s + Ip^k) + ^j^{M^, + (M^ - s + 2p2kf 



+ |c6p(M4 - s + 2p2fc)2 + \c7W2kp1] 



+ ^K[c5C^m2^(M4 - s + 2p2fc) + C54A4(2fcpi) + C6C^(M4 - S + 2p2k){2kpi)] + ^[|C5|"M^^ 



(19) 

We have here used a combination of the notation from Ref . [2|l25p26j so that the definitions 
which follow should hold true: 

®i = {q+ Pj)'^ - nij + ie = + 2qpj + pj - rrij + ie = + 2qpj + r^, 

Api,j =Pi-pj, 

DqIJ) = 4-pomt scalar function obtained from 5-point scalar function by omitting denominator Tij, 
Co{i,j) = 3-point scalar function obtained from 5-point scalar function by omitting denominators 
and D j , i ^ j, 

where we also follow the Passarino-Veltman [2] notation of the 4-point one-tensor integral, 
Dfj,{j), obtained by omitting denominator Dj from the corresponding 5-point one-tensor 
integral with 

, where the 4- vectors {p{j)i} are then determined in accordance with Ref. [2], with the 
understanding that p(j)4 is only non-zero if it is necessary to shift the q integration variable 
by it to reach the standard form of the respective Passarino-Veltman representation. This 
expression for Eq does not have problems with Gram determinant type denominators. 

To further exhibit the magic in the polarization vector spinor representation under 
display here, we record as well the results for Fig. (la) and (lb) that one needs to add to 
our result for Fig. (Ic) to get a gauge invariant result: 

= 0, by 'magic' 

<-^-^ = ^'-^ ^(Pi +P^-P^-P^-'^ Kk-\p,+ ><k- \p2+ > fep 
where the numerator A^' is given by 

N' = (< P2 - \Pl+ > ai+ < P2 - \P2+ > bi) (< pi - \p2+ >< P2 + \p'l- > 

- <pi-\k+ >< k + \p[- >) + (< P2 - \pi+ > ai+ < P2 - \p2+ > bi) (21) 
[{-2pi(j)2 - k))a4 + 02 < Pi - \k+ ><k + \p2- >] 
with the definitions 

oi = (1 + ai){2pip[) + a3<p2 + \p'i- >< p'l - \pi+ > 

bi = a2 < P2 + \p'i- >< p'l -\pi+> +a4{2pip'i) ^^^^ 
ai =< pi - \P2+ > [(1 + ai) <Pi + \p'i- > +a3 < P2 + \p'i- >] 
bi =< pi - \P2+ > [a2 <P2 + \p'i- > +OLi <pi + \p'i- >]. 
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Again, this gives immediate reduction to the known scalar functions with considerable 
reduction in the number of terms requiring evaluation compared to the usual trace over 
fermion lines method when one computes the respective contribution to 2^M.*^Ai^^. The 
complete phenomenology of our results for the process in Fig. [T]will appear elsewhere [36] . 

It is important to explain the difference between what we have done here and what was 
done in Refs. [3l|4l|40l|41] . We do this in turn in a somewhat reverse chronological order. 
In Ref. [IQ], the representation of the loop variable in a basis of light-like 4- vectors is used 
to construct a recursion relation between one-loop n-point tensor integrals of differing 
rank whereas in Ref. [H] the spinor representation of the external tensor coefficient of 
a massless n-point tensor one-loop integral is used to reduce the rank of that integral 
iteratively to allow numerical implementation, using Dirac matrix methods. In both cases, 
the square roots of the Gram determinants appear in the denominators of the resulting 
representations. In our approach, explicit kinematics allows direct Chinese magic action 
in the complete amplitude contribution's evaluation directly to the lower point functions 
without Gram determinant factors to be computed in our denominators. No iteration is 
necessary and Chinese magic reduces considerably the number of terms in our final result. 
Such action is not present in Refs. [IHlllI]. In Refs. the representation of the n-point 
amplitude at one-loop starts from its integrand iV(g)/(S?o ■ ■ ■ ^n-i) with an expansion of 
the numerator N{q) in powers of the denominators with coefficients that split into a 
part that is independent of q and a part that integrates to zero with the understanding that 
the integration measure is in general in d dimensions whereas the function N{q) is defined 
for q in 4-dimensions. We refer to this representation as the OPP representation after the 
authors in the first paper in Refs. Various methods for adding in the so-called missing 
rational terms generated by the mismatch between the 4-dimensional q in N and the d- 
dimensional q in the {Sj} are given in Refs. including the generalized d-dimensional 
unitarity that treats the full d-dimensional unitarity realization of the OPP representation. 
In all of these works, N[q) or N{q) is treated as a given and no procedure for exploiting 
Chinese magic to simply it at the loop momentum level is considered. Moreover, the need 
to add in rational terms is an essential part of the procedure, whereas, as we see in our 
result (fT3l) . we do not have such an issue in our approach - we get the complete answer 
with methods that operate entirely in 4-dimensionqj. More importantly, inverse powers 
of Gram-type determinants appear in the coefficients in the representation of N so that 
issues of numerical stability obtain, whereas as we show above our approach does not lead 
to such factors so that it should be more stable. Finally, the procedure for determining 
the coefficients in the representation of involves solving the algebraic problem for the q 
values at which 4 , then 3, then 2, and finally 1 of the vanish(es). This means that, 
in general, complex values of q are required and this forces the evaluation of N{q) at such 
unphysical 4-momenta. Our approach avoids this issue altogether as we carry our entire 
calculation out in the 4-dimensional real virtual loop momentum space. We then provide 

^If one wants to apply our method to lower point amplitudes that are UV divergent, in renormalizable 
theories one should use the known counter-terms for those divergences to render the amplitudes finite 
first and then apply our 4-dimensional methods to the UV finite subtracted amplitudes. 
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a completely physical cross check on the methods in Refs. Similarly, the approach 
in Refs. [3] also takes the integrand as a given and constructs the respective amplitude 
from unitarity-based on-shell (recursion) relations, where the authors in Refs. [3] are 
able to get both the cut-constructable and the rational parts of the amplitudes with 
such methods. Again, there is no exploitation of Chinese magic to simply the amplitude 
at the loop variable level, the amplitude construction uses 4-particle cuts that have in 
general complex 4-momenta as their solutions so that wave functions are evaluated at 
such unphysical momenta, and the solution of these on-shell relations generally introduces 
troublesome kinematic factors in the denominators of the representation so that numerical 
stability cannot be assured. Our approach avoids all of these problems and affords again 
a completely physical cross check on this approach as well. 

The complete analytical result for the amplitude in Fig. [1] will be presented else- 
where [36] • Here, we have shown that the use of Chinese magic in the virtual loop 
momentum can reduce considerably the amount of algebra required for stable, efficient, 
manifestly physical computation of higher point virtual corrections with general mass 
scales, as they are needed for exact amplitude-based resummed MC event generator real- 
ization. 
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